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a 
Abstract 
Produced water re-injection (PWRI) is often the safest and most economical method for disposal of produced water in the 
oil industry. Two key issues that affect the management of PWRI are the formation damage and the constrained pumping 
pressure at the wellhead. A simulator was developed to handle the design of single-zone or multi-zone water injection in 
multi-layered reservoirs. The simulator can accommodate both vertical and horizontal wells operated under matrix 
and/or fractured regimes. It is also able to account for the impact of formation damage and user-defined wellhead 
pressure constraints. Results obtained from the simulator showed good agreement with known injection behaviors. For 
vertical wells, injection conformance depends on KH (permeability-thickness) and the minimum horizontal stress; in the 
case of multi-fractured horizontal wells, the outermost fractures (those near the tip and the heel of the horizontal well) 
are longer than the fractures in the middle. Lastly, by constraining the maximum allowable surface pressure, frictional 
pressure drops in both the wellbore and fracture cause the injection rate to decline, which in turn affects both the 


fracture geometry and the maximum disposal volumes. 


Keywords: Produced water re-injection; Formation damage; Hydraulic fracture; Fracture simulator; Injectivity decline; 


Constraint wellhead pressure 
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Introduction key factors, including water quality, formation 
compatibility, reservoir properties, and _near-well 

Produced water re-injection (PWRI) provides an geomechanics [1-5]. Moreover, many PWRI wells 
environmentally acceptable solution to the waste water assumed to be operating under matrix injection 
disposal, reservoir pressure maintenance and enhanced conditions have actually become hydraulically fractured. 
oil recovery. PWRI performance depends on a number of Our results indicate that this transition from matrix to 
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fractured injection is often essential to the long-term 
maintenance and successful operation of the PWRI well. 


PWRI performance depends critically on evolution of 
near-well formation damage. Both chemical and 
mechanical processes may lead to a loss of matrix 
injectivity even in highly permeable formations. Two 
mechanisms are often to blame: particulate matter 
carried within the injected water and _ chemical 
precipitates formed within the wellbore or reservoir 
during injection may become deposited on rock surfaces 
to form a heterogeneous and highly compressible filter 
cake that clogs the pore throats and/or fracture surfaces 
[6-12]. In each of these scenarios, the permeability near 
the well decreases and injectivity begins to decline [13]. 


During the period of formation damage, surface 
pressures are observed to increase as the desired flow 
rate is maintained. (Note that the maximum surface 
pressure may act as an important constraint.) If the 
bottom-hole pressure rises above the breakdown 
pressure, a fracture is created and the well injectivity 
subsequently increases: the newly formed fracture allows 
for a rapid increase in the accommodated fluid volume, 
and the fracture’s surfaces enable the injected water to 
leak off into undamaged rock further from the well. 
Continued PWRI operation causes the fracture to grow 
and well performance to improve. Indeed, advanced step- 
rate and fall-off tests are able to determine this transition 
from matrix to fractured injection, as well as to determine 
the influence of water properties, leak-off and reservoir 
permeability on each mode of injection. 


@FRAC2D Simulator 


The goal of this section is to outline a newly developed 

PWRI simulator that satisfies these requirements: 

e Account for the damage that results from solids and 
oil in the injected water 

e Account for the possibility of a constrained well-head 
pressure and allow injection under constant flow rate 
and under constant surface pressure 

e Permit matrix and/or fractured injection into multi- 
layered formations for both vertical and horizontal 
wells 

e Include effects of “stress shadowing” on _ the 
propagation of multiple fractures in horizontal wells 

e Incorporate “flow partitioning” within the wellbore to 
accurately account for dynamic changes in fluid 
pressures and rates in the case of multiple 
propagating fractures 


Simulator Background 


Numerical models of matrix injection and hydraulic 
fracturing under conditions of variable formation damage 
are described in what follows. Models of both deep-bed 
filtration (which consist of equations for mass balance of 
suspended and deposited particles) and kinetics-based 
particle capture by the rock are also needed to adequately 
predict how the permeability will be modified during 
injection [14-16]. Techniques for implementing surface 
pressure constraints inside the flow simulator are also 


discussed. 


Table 1 indicates how each of the model input 
parameters is related to values that can be recovered 
from well tests, geologic data, well logs and core analyses. 
Simulator results include both time series and fracture 


geometries. 


Input 


Output 


Well Data 


Wellbore radius 


Tubing length, 
roughness, ID 


Injection Parameters 


Injection fluid 
temperature 


- Injection schedule (rate, 
time) 


- Injection fluid rheology 


Completion Data 


- Perforation top and 
bottom 


Zone deviation, 
thickness 


Flow rate distribution 
between the layers and 
between multiple 
perforations 


Wellhead pressure 
(WHP) 


Bottom Hole Pressure 
(BHP) 


Mechanical Properties 


Minimum horizontal 
stress 


Young’s modulus 


Poisson’s ratio 


Fracture toughness 


Rock surface energy 


Injectivity Index (II) 


Reservoir Properties 


Reservoir pressure 


- Reservoir temperature 


Porosity and 
permeability 


Flow contributions for 
both matrix and fracture 
injection 
Fracture length with 
consideration of poro- 
elastic and thermo-elastic 
effects 
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Table 1: Input and Output Parameters for the Simulator. 
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Transition from Matrix to Fracture Injection 


Prior to fracture initiation, the simulator uses multi- 
zone KH (permeability-thickness) analysis, wellbore 
friction, and reservoir pressure to predict injectivity 
conformance profiles along the wellbore. The transition 
from matrix to fractured injection is controlled by the 
following condition [17] 


Pup 20,+ a (1) 
2(1—v")r,, 


Where 0, is total stress perpendicular to the fracture: 


O, =(O x) pin + AG yp + Ao", (2) 
Here: 
Piwt = Flowing bottom-hole pressure, psi 

V = Poisson’s ratio 

U = Rock surface energy, ft-lbf/sq in. 

E = Young's modulus, psi 

(c,, Ve = Minimum in-situ horizontal stress, psi 
Tw — Radius of the wellbore, ft 


The simulator models each fracture as if it were 
embedded within an elliptical (multi-layered) cylindrical 
reservoir. Each layer is assumed homogeneous. Multiple 
concentric reservoirs can be included to determine the 
extent of formation damage far from the fracture. In this 
case, the inner cylinder contains the fracture and is 
allowed to become damaged; a cylindrical reservoir 
containing the inner cylinder captures the effect of the 
background (unaltered) reservoir outside the extent of 
the inner zone. Dynamic equilibration between the 
fracture’s internal pressure and far-field pressure leads to 
additional terms in the total stress that acts on the 
fracture’s evolution: 


Temporal change (final - initial) in average 
interior stress perpendicular to the major axis of 
Ao", = the ellipse resulting from a pressure difference (p 
- Pg) between the elliptical cylinder and the 
surroundings 
Temporal change (final - initial) in average 
interior stress perpendicular to the major axis of 
the ellipse resulting from a temperature 
difference (T - Tp) between the elliptical cylinder 
and the surroundings 


AG = 
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Fracture Propagation Model 


Once the fracture has been initiated, the pressure at the 
fracture tip is determined using: 


MUE 


—,— 3 
2(1-v*)r, ” 


pi =o0,+ 


The difference between the wellbore pressure and the 
average pressure at the face of the fracture is determined 
using an expression that relates the fracture geometry, 
uniform leak-off and internal flow (assumed to be 
laminar) [18]: 


P, = Pup — Ap f (4) 
where 
3 4 
Ap Ly HG : 
4 
and 
ly = Water injection rate, BBL/day 
hos Fluid viscosity, cp 
lf = Fracture half length, ft 
h 7 Reservoir thickness, ft 
G = Elastic shear modulus, psi 
Radius of the extending edge of the fracture 
rf = (equals the smaller of the two values: Lyor h/2), 


ft 


The fracture width W is controlled by the fracture net 
pressure, Pnet: 


W = F(Pret) (6) 


Expressions for the nonlinear function F are 
numerically complex, but analytical examples can be 
found in a variety of sources. 


Plugging Damage Model 


Formation damage considered in the simulator focuses 
on the reduction of permeability resulting from injection 
of solids. Three pore-scale scenarios are explored in 
Figure 1: 
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Cen 
C 


External filter cake 


Internal filter cake Non-bridging solids 


ccrpies bs soda > 33% dpore 33% > dicta > 14% dyore —_duodd < 14% pore 


Low velocity (<2 
cm/min interstitial) 


sod > 33% dpore 33% > dsoila > 7% dpore soitd < 7% dpore 


Figure 1: Damage mechanism due to solid injection [19]. 


1. Plugging of pores by solids when pore throat range is 
3 times the particle size or less. The result is a 
reduction in permeability that can be reflected as a 
reduction in porosity from the porosity-permeability 
correlations. 

2. Settling of solids into the pores when the pore throat 
range is three to seven times the particle size. The 
result is a reduction in porosity, which is reflected as 
a reduction in permeability from the porosity- 
permeability correlations 

3. Solid particles that pass through the pores with a 
pore throat range at least seven times the particle 
size. In this case, the solid particles pass through the 
pore throat with no possibility of plugging and no 
reduction in permeability is expected. 


The pressure increase due to skin damage around the 
wellbore is calculated following [20]. 


AP. = 1, Ll, 4, (7) 
where 
R, =R,, +R, —|In(R, /7,)/ (22hK)| (8) 


and Rincis the internal filter cake obtained from 


= | (9) 


The external filter cake, R,, is obtained from: 
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_ ints, Ar, —2.)) 


R. (10) 
27nK . 
Where 
Ra 7 The damaged radius, ft 
K(r) = The permeability at each radius, mD 

K _ The permeability of the external filter cake, 
oc mD 

he 7 Filter cake thickness, ft 


Porosity Reduction Model for Oil-in-Water 
(OIrw) 
In this model, the oil in water is converted to an 


equivalent solid concentration and assumed to behave as 
if it were additional solid particulate: 


C= OIW* 0.14/2.1 (11) 
Here 
C The equivalent solid concentration of oil in 
~ water, ppm 
OIW = Oil in water concentration, ppm 


The empirical correlation in Equation (11) is one example 
of a functional form that was developed by Advantek 
International as part of comparisons and best-fit analyses 
to field data using @FRAC2D from multiple PWRI 
projects. Other examples have also been used [21-22]. 


@FRAC2D Validation 


Numerous examples have been used within Advantek 
International to validate the simulator. Three particular 
results will be discussed here: 

Case #1: Comparison to published results 

Case #2: Vertical well with thief zones 

Case #3: Horizontal well with stress shadowing 


Case #1: Comparison to Published Results 


Table 2 summarizes the input data for the injected well 
and the geomechanical parameters for different layers for 
a test case that had previously been published in the 
literature [23]. The injection rate for the well was 25,000 
bpd. The solid concentration was 5 ppm with an average 
diameter of 5 microns and the particle density was 2.3 
gm/cm3. 
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Layer Layer Layer 
L1 L2 L3 
Depth (mid Layer), ft | 16925 | 16932.5 | 17027 
peer 180 180 180 
Temperature, F 
Thickness, ft 30 35 45 
Porosity 0.27 0.28 0.3 
Permeability, md 407 529 687 
Zoune Modus, 0.165 0.12 0.12 
Mpsi 
Poisson’s Ratio 0.25 0.28 0.28 
Min, Marontal 9500 | 10200 | 10200 
Stress, psi 


Table 2: Input parameters for the published case. 


Injection Time (Da' 


Figure 2A: Fracture growth obtained from the 
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Figure 2B: Fracture growth obtained using @FRAC2D. 
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Figure 3A: Fluid distribution obtained from the 
published case. 


Comparisons of injection conformance were also made. 
Figure 3A, which contains the published results, indicates 
that approximately 90% of the injected fluid flowed into 
Layer #1. 


published case. 


Comparisons between the published results and the 
results from @FRACZ2D are very good. Figure 2A, which 
contains the published results, indicates that the fracture 
grew within Layer #2 during the first day of injection and 
that the fracture length extended to approximately 46 ft. 
Tiny fractures were observed in the other layers. Figure 
2B, which contains the @FRAC2D results, indicates 
approximately the same results: the fracture length grew 
to approximately 47 ft and exhibited the same growth 
behavior. 
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Figure 3B: Fluid distribution obtained using 
@FRAC2D. 
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Figure 3B, which contains the @FRAC2D results, 
indicates essentially the same qualitative and quantitative 
behaviors. 


The simulator input and formation parameters for the 
vase of a vertical well are summarized in Tables 3 and 
Table 4. Three distinct scenarios are described below. 


Case #2: Vertical Well 
Injection Time 25 years 
Injection Rate 35,000 bpd 
Tubing Roughness 0.0001 
Tubing Head Pressure 3,500 psi 
Minimum Horizontal Stress, Shale omin 0.8 psi/ft 
Total Solid Concentration, TSS 5 ppm 
AT -60 °F 
Oil in Water, OIW 10 ppm 
Poisson Ratio 0.24 
Young’s Modulus for Shale, Eshate 55,0000 psi 
Young’s Modulus for Sand, Esana 85,0000 psi 
Sand Linear Coefficient of Thermal Expansion 6.5E 
Shale Linear Coefficient of Thermal Expansion 4.5E6 
Tubing ID 6 inch 
Table 3: The input parameters for the @FRAC2D. 
Formation Layer | Layer | Layer | Layer wee Ree ee eee 
Parameters 1 2 3 4 16000 
K, md 2000 1000 4000 8000 a 
h, ft 15 60 40 10 Z 10000 
Omin, psi/ft 0.6 0.55 0.58 0.6 i 7 
Top, ft 9415 9510 9620 9700 : 4600 
Bottom, ft 9430 570 9660 9710 zone 
Pressure, psi/ft} 0.445 0.445 0.445 0.445 *0 1000 2000 a ee ee 


Table 4: The formation parameters. 


Matrix Injection with No Damage: In this scenario, it 
was assumed that no damage would occur during 
injection (i.e., 100 % of the total injection solid will pass 
through the formation.). The wellhead pressure was 
assumed to be 3,500 psi. 
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Figure 4: Fluid distribution under matrix injection. 


Figure 4 shows that the injection rate was distributed 
between four layers based on the various KH 
(permeability-thickness) values, as expected. Layer #3, 
which had the highest KH value, received the highest flow 
rate, while Layer #1 received the lowest flow rate. The 
injection rate was constant for the four layers during the 
injection time because there was no damage. This 
scenario confirmed that when no damage or fracturing 
occurs, flow portioning between the layers nearly mirrors 
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the KH distribution. Differences from this distribution are 
due to the effects of wellbore friction and differing 
reservoir pressures, which slightly decreases the flow to 
the deeper layers, as shown in Table 5. 


Layer Layer Layer Layer 
1 2 3 4 
0, 
ane 9.1 18.2 48.5 24.2 
total 
Q % of total 9.3 19.1 48.2 35.5 


Table 5: Fluid distribution is a mirror of the KH in case of 
no damage or fracture. 


Distributed Damage: In this scenario, it was assumed 
that 80% of solids pass through Layer #1, 60% of solids 
pass through Layer #2, 90% of solids pass through Layer 
#3, and 100% of solids pass through Layer #4. The 
wellhead pressure was again assumed to be 3,500 psi. As 
damage accumulates, the undamaged layer acts as a thief 
zone which accumulates all of the flow, as shown in 
Figure 5. 
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Figure 5: Undamaged zone acts as a thief zone. 
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Figure 7: Fluid distribution and formation plugging. 


If it was assumed that 80% of the injected solids pass 
through the four layers and the wellhead pressure was 
3500 psi. Damage accumulated in each of the layers. Layer 
#2 fractured after 430 days, but the fracture propagated 
only slightly (until it reached 2 ft long), as shown in 
Figure 6, Figure 7 indicates that as the fracture 
propagated, most of the fluid flowed in this layer (90% of 
the injection fluid). 
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After 2,705 injection days, the fracture in Layer # 2 
became completely plugged and as a result, the injected 
fluid distributed again through the other layers, which 
caused each to hydraulically fracture. As the damage 
accumulated in the shortest fracture, Layers #1 and #4 
were plugged after 4,030 days. This resulted in all of the 
fluid flow redistributing to Layer #3. After 4,670 days, 
Layer #3 became completely plugged and injection was 
stopped due to complete formation plugging. 


Damage with Constrained Surface Pressure: In this 
case, it was assumed that 80% of the injected solids could 
pass through the formation and that the wellhead 
pressure was assumed to be only 1,350 psi. 
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Figure 8: Impacts of imposed constraints on the 
surface pressure. 


Figure 8 indicates that after 110 days, the required 
well-head pressure reached the maximum value (of 1,350 
psi) because of damage - as a result, the injected flow rate 
declined and no fracture were occurred in any layers. 
Oscillations in the predicted wellhead pressure arise 
because of the particular method that was used to enforce 
the constraint; new penalty methods are being developed 
to ensure that the constraint is met exactly at each time 
step, which results in a smooth pressure response. 


Case #3: Horizontal Well 


The term “stress shadow” corresponds to the stress 
field generated by a hydraulic fracture. This stress field 
impacts neighboring fractures, faults, the surrounding 
reservoir (and pore pressure field), as well as the 
wellbore itself. The increase in stress is greatest at the 
fracture face and fracture tip, but the stress variation 
extends out into the reservoir for hundreds of feet. The 
impact of stress shadows is essential for optimizing frac 
job planning in horizontal wells, since fractures are 
initiated simultaneously at multiple perforation clusters. 

@FRACZ2D is able to determine how the stress field 
generated by each fracture (including the roles of leak-off 
and net pressure) will impact each of the neighboring 
fractures’ geometries. The results from a simple example 
of injection into a horizontal well are shown in Figure 9. 
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Figure 9: Effect of the stress shadow on the fracture 
length for a horizontal well. 
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In this case, the fractures were initiated uniformly 
along the wellbore. The longest fractures were found 
close to the heel and the toe of the well (at 72 and 75 ft, 
respectively), while the shortest fracture was found in the 
middle (52 ft in length). In general, increased compressive 
stresses from a dominant fracture were found to prevent 
the initiation of nearby parallel fractures, which provided 
a natural diversion mechanism along the wellbore. If 
fractures were initiated too close together, stress 
shadows tended to decrease fracture growth along the 
mid-section of horizontal wellbores and to encourage 
fracture growth at the heel and toe of the wellbores 
(similar to results found elsewhere [24].) 


Conclusion 


Results from @FRACZ2D have been shown to effectively 
and accurately re-produce expected hydraulic fracture 
behaviors for a range of cases. The impact of both matrix 
and fracture injection can be understood, as can the role 
of thief zones and formation damage on_ injection 
conformance. @FRACZ2D also captures formation damage 
caused by solids and solid settling, and oil in water during 
both matrix and fractured injection. Simulation results 
also compare very well to results from other published, 
benchmark cases. 


@FRAC2D, which considers stress changes due to 
thermal and poro-elastic effects, has also been shown to 
predict fracture containment (vertical growth) and stress 
shadowing (for horizontal wells with multiple, 
propagating fractures). The simulator can _ also 
accommodate various penalty methods to enforce facility 
constraints such as a maximum allowable surface 
pressure and constant rate or constant pressure. In 
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conclusion, @FRACZD is a fast, numerical tool that allows 
an operator to understand the impact of injection fluid 
properties on injectivity, pump horsepower 
requirements, surface treatment needs, maximum 
disposal rates and volumes. Advantek International 
continues to improve the simulator’s robustness and 
underlying models as it is used for PWRI and other field 
management projects. 
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